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Abstract 

Let (A(n))^gj!^ be a sequence of i.i.d. topical (i.e. isotone and addi- 
tively homogeneous) operators. Let x{n,xo) be defined by x(0,xo) = 
xo and x{n, xq) = A{n)x{n — 1, xq). This can modelize a wide range of 
systems including, task graphs, train networks, job-shop, timed digital 
circuits or parallel processing systems. 

When (A(n))„gj^ has the memory loss property, we use the spectral 
gap method to prove limit theorems for x(n, xq). Roughly speaking, 
we show that x(n, xq) behaves like a sum of i.i.d. real variables. Pre¬ 
cisely, we show that with suitable additional conditions, it satisfies 
a central limit theorem with rate, a local limit theorem, a renewal 
theorem and a large deviations principle, and we give an algebraic 
condition to ensure the positivity of the variance in the CLT. When 
A(n)s are defined by matrices in the (max, -t-)semi-ring, we give more 
effective statements and show that the additional conditions and the 
positivity of the variance in the CLT are generic. 


Introduction 

An operator A : is called additively homogeneous if it satisfies 

A{x + al) = A{x) + al for all x G and a G M, where 1 is the vector 
(1, • • •, 1)' in It is called isotone ii x < y implies A{x) < A{y), where 
the order is the product order on M'^. It is called topical if it is isotone and 
homogeneous. The set of topical operators on will be denoted by Topd- 
We recall that the action of matrices with entries in Mmax = M U {—oo} 
on is defined by {Ax)i = maXj(Ay -f- Xj). When matrix A has no line 
of —oo, the restriction of this action to defines a topical operator, also 
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denoted by A. Such operators are called (max, +)operators and composition 
of operators corresponds to the product of matrices in the (max, +)semi-ring. 


Let (yl(n))^gj^ be a sequence of random topical operators on 
x{n,xo) be dehned by 

a:(0,Xo) = xo 

x{n,Xo) = A{n)x{n — l,Xo). 


Let 


( 1 ) 


This class of system can modelize a wide range of situations. A review of 
applications can be found in the last section of |BM98j . When the x{n,.) are 
daters, the isotonicity assumption expresses the causality principle, whereas 
the additive monotonicity expresses the possibility to change the origin of 
time. (See J. Gunawardena and M. Keane |GK95j . where topical applica¬ 
tions have been introduced). Among other examples the (max,-|-)case has 
been applied to modelize queuing networks (J. Mairesse |Mai97] . B. Hei- 
dergott |Hei00j ). train networks (B. Heidergott and R. De Vries |HdV01] . 

H. Braker |Bra93j i or Job-Shop (G. Gohen and ah |GDQV85l ). It also com¬ 
putes the daters of some task resources models (S. Gaubert and J. Mairesse |GM98j i 
and timed Petri Nets including Events graphs (F. Baccelli [Bac92j ) and 1- 
bounded Petri Nets (S. Gaubert and J. Mairesse |GM99] ). The role of the 
max operation is synchronizing different events. For devlopements on the 
max-plus modelizing power, see F. Baccelli and ah BGOQ9^ or B. Heider¬ 
gott, G. J. Olsder, and J. van der Woude |HOvdW06] . 

We are interested in the asymptotic behavior of x{n,.). It follows from 
theorem 11.11 that - maxj a:j(?T., Vq) converges to a limit 7 . 

In many cases, if the system is closed, then every coordinate Xi{n,Xo) 
also converges to 7 . The value 7 , which is often called cycle time, is the 
inverse of the throughput (resp. output) of the modelized network (resp. 
production system), therefore there has been many attempt to estimate it. 

(J.E. Gohen [Goh 88 j . B. Gaujal and A. Jean-Marie |GJM98] . J. Resing and 
al. |RdVH^9n] ) Even when the A(n)’s are i.i.d. and take only hnitely many 
values, approximating 7 is NP-hard (V. Blondel and al. |BGT00j i. D. Hong 
and its coauthors have obtained f |BH00aj . |BH00b] . [GHOOj ) analyticity of 
7 as a function of the law of A(l). In this paper, we prove another type of 
stability, under the same assumptions. 

We show that under suitable additional conditions, x(n,.) satishes a cen¬ 
tral limit theorem, a local limit theorem, a renewal theorem and a large 
deviations principle. When the A{n) are (max, -|-)operators we give more 
explicit results. Those results justify the approximation of 7 by ^Xi(n,Xo) 
and lead the path to conhdence intervals. 
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Products of random matrices in the usual sense have been intensively in¬ 
vestigated. Let us cite H. Furstenberg |Fur63j . Y. Guivarc’h and A. Raugi |GR85j 
or I. Ya. Gol'dsheid and G. A. Margulis |GdM89] . The interested reader 
can hnd a presentation of this theory in the book by Ph. Bougerol and 
J. Lacroix |BL85] . We investigate analogous problems to those studied by E. 

Le Page [LP82j . but for matrices in the (max, -|-)semi-ring and more generally 
for iterated topical operators. 

This article is divided into three parts. First we present the model of 
iterated topical operators, including a short review of known limit theorems 
and a sketch of the proof of our results. Second we state our theorems and 
comment on them. Finally we prove them. 

1 Iterated topical operators 

1.1 Memory loss property 

Dealing with homogeneous operators it is natural to introduce the quotient 
space of by the equivalence relation ~ dehned by x ~ ?/ if x — |/ is 
proportional to 1. This space will be called projective space and denoted by 
^®max- Moreover x will be the equivalence class of x. 

The application x h-^ {xi — Xj)i<^j embeds onto a subspace of R ^ 2 ' 

^ _ d{d-l) 

with dimension d—l. The inhnity norm of R 2 therefore induces a distance 
on PRmax which will be denoted by 6. A direct computation shows that 
S(x, y) = maxi(a:i — yi) + ma.Xi{yi — Xi). By a slight abuse, we will also write 
5{x, y) for 6(x,y). The projective norm of x will be \x\'p = 5{x, 0). 

Let us recall two well known facts about topical operators. First a topical 
operator is non-expanding with respect to the inhnity norm. Second the 
operator it dehnes from PR^ax to itself is non-expanding for 5. 

The key property for our proofs is the following: 

Definition 1.1 (MLP). 

1. A topical operator A is said to have rank 1, if it dehnes a constant 
operator on PR^ax • does not depend on x G R'^. 

2. The sequence (A(n))^gj^ of Top^-valued random variables is said to 
have the memory loss (MLP) property if there exists an N such that 
A{N) ■ ■ ■ A(l) has rank 1 with positive probability. 

This notion has been introduced by J. Mairesse [Mai97j . the A{n) being 
(max, -|-)operators. The denomination rank 1 is natural for (max, -|-)operators. 
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We proved in [Mer04j that this property is generic for i.i.d. (max, +)operators: 
it is fnlhlled when the snpport of the law of A(l) is not inclnded the union 
of hnitely many affine hyperplanes. 

Although this result could suggest the opposite, the MLP depends on the 
law of A(l), and not only on its support : if is an i.i.d. sequence 

with the support of U{1) equal to [0,1], and A{n) are the (max, +)operators 
dehned by the matrices 


A{n) 


-U{n) 0 \ 

0 -U{n) ) ’ 


then (A(? 7 ,))^gj^ has the MLP property iff F{U{n) = 0) > 0. 

The weaker condition that there is an operator with rank 1 in the closed 
semigroup generated by the support of the law of A(l) has been investigated 
by J. Mairesse for (max, +)operators. It ensures the weak convergence of 
x(n,.) but does not seem appropriate for our construction. 


1.2 Known results 

Before describing our analysis, we give a brief review of published limit the¬ 
orems about x{n,X^). 

There has been many papers about the law of large numbers for products 
of random (max,-|-)matrices since it was introduced by J.E Cohen |Coh 88 ] . 
Let us cite F. Baccelli |Bac92j . the last one by T. Bousch and J. Mairesse [BMObj 
and our PhD thesis [MerOb] (in French). The last article proves results for a 
larger class of topical operators, called uniformly topical. 

J.M. Vincent has proved a law of large number for topical operators, that 
will be enough in our case : 

Theorem 1.1 ( |Vin97] ). Let be a stationary ergodic sequence of 

topical operators and an W^-valued random variable. If A{1).0 and X^ 
are integrable, then there exists 7 and 7 m M such that 

maXixAn, X^) _ 

hm- = -y a. s. 

n n 

mhiiXi{n,X^) 

hm- = 7 a.s. 

n n — 

F. Baccelli and J. Mairesse give a condition to ensure 7 = 7 , hence the 

r x{n,X^) 

convergence oi 

Theorem 1.2 ( |BM98j ). Let (A(n))^gj^ be a stationary ergodic sequence of 
topical operators and X^ an W^-valued random variable such that A(1).0 and 
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are integrable. If there exists an N, such that A{N) ■ ■ ■ A(l) has a bounded 
projective image with positive probability, then there exists 7 m R such that 


lim 


x{n, 


n n 


'jl a.s. 


In this case 7 is called the Lyapunov exponent of the sequence. We notice 
that the MLP property implies a bounded projective image with positive 
probability. 

The following result has been proved by J. Mairesse when the A{n) are 
(max, +)operators, but can be extended to topical operators with the same 
proof. It will be the key point to ensure the spectral gap. 

Theorem 1.3 (Mairesse |Mai97j L If the stationary and ergodic seguence 
(y4(n))^gg of random variables with values in Topd has memory loss property, 
then there exists a random variable Y with values in PR^ax '■ = 

A{n) ■ ■ ■ A{\)Y is stationary. Moreover 


lim P (dxo, Yn 7 ^ x{n, xq )) = 0. 

n—^oo 

In particular x{n, Xq) converges in total variation uniformly in xq. 

The law of Y is called the invariant probability measure. 

To end this section we mention two limit theorems, which are close to 
ours, but obtained by different ways. We will compare those results to ours 
in section 12.31 

With a martingale method J. Resing and ah |RdVH'*~90] have obtained a 
central limit theorem for x{n, when the Markov chain x{n, .) is aperiodic 
and uniformly <I)-recurrent. The theorem has been stated for (max, +)operators, 
but it should make no difference to use topical ones. 

With a subadditivity method, F. Toomey |Too02j has proved a large 
deviation principle for x{n,Xo) when the projective image of A{N) ■ ■ -^(1) 
is bounded. 


1.3 Principle of the analysis 

From now on, (y 4 (? 7 ,))^gj^ is an i.i.d sequence of topical operators with the 
MLP property. 

The first step of the proof is to split our Markov chain x{n ,.) into another 
Markov chain and a sum of cocycles over this chain, following what E. Le 
Page made for products of random matrices. For any topical function 0 from 
R"^ to R, (j){Ax) — 4>{x) only depends on A and x. Therefore (j){x{n ,.)) — 
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(j){x{n — 1,.)) only depends on A{n) and x{n — 1, .)• Since can be seen 

as an hyperplane of x{n,.) can be replaced by {(j){x{n, .)),x{n,.)). (cf. 
lemma 13.31) 

According to theorem 11.31 we know that x{n ,.) converges. On the other 
hand, by theorem 11.21 xin. X°). goes to inhnity (if 7 7 ^ 0 ) in the direction 
of 1, so (f>{x{n ,.)) ~ yn. We investigate the oscillations of (f>{x{n ,.)) — yn. 
Interesting 0’s are dehned by 0(x) = Xj, 0(a:) = max^Xj, 0(x) = minjXj. 

The second step is to prove the spectral gap for the operator dehning the 
Markov chain {A{n),x{n — 1, and apply the results of [HHOl] et [HHOdj 

that give limit theorems for (j){x{n, X^)) — 0(W°) — yn. The spectral gap 
follows from the convergence of x{n ,.), just like by E. Le Page |LP82j . 

We use two series of results. The hrst series are taken from the book 
by H. Hennion and L. Herve [HHOl] that sums up the classical spectral gap 
method developed since Nagaev |Nag57| in a general framework. To apply 
it we demand integrability conditions on sup^, | 0 (A(l)x) — 0 (x)| to have a 
Doeblin operator on the space of bounded functions. The second series are 
taken from the article [HH04j that is a new rehnement of the method in 
the more precise framework of iterated Lipschitz operators. Since our model 
enters this framework, we get the same results with integrability conditions 
on A(1)0 that ensures that the Markov operator satisfy a Doeblin-Fortet 
condition on functions spaces dehned by weights. The comparison between 
the two series of results will be made in section 12.31 

2 Statement of the limit theorems 

2.1 General case 

From now on, we state the results that we will prove in section [31 

For local limit theorem and for renewal theorem we need non arithmetic- 
ity conditions. There are three kind of non arithmeticity, depending if the 
theorem follows from [HHOl] or [HHOl] . We will denote them respectively by 
(weak-) non arithmeticity and algebraic non arithmeticity. When d = 1 they 
fall down to the usual non arithmeticity condition for real i.i.d. variables. 
Algebraic non arithmeticity will be dehned before the statement of LLT, but 
other non arithmeticity conditions will be dehned in section [3] once we have 
given the dehnitions of the operator associated to the Markov chain. Unlike 
algebraic non arithmeticity, they depend on the 2 -uple ((A(n))^gj^, 0 ), which 
will be called ’’the system”. 

Let (A(n))^gj^ be an i.i.d sequence of topical operators with the MLP 
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property. The sequence (a;(n, is defined by equation ([T]) and 7 is the 

Lyapunov exponent defined by theorem II.2[ 

Since the topology of the uniform convergence over compact subset on 
Topd has an enumerable basis of open sets, the support of measures on it is 
well defined. We denote by Sa the support of the law of ^(1) and by Ta the 
semi-group generated by Sa in Topd- 

Theorem 2.1 (CLT). Let (y 4 (? 7 ,))^gp^ be an i.i.d sequence of topical operators 
with the MLP property and an M.'^-valued random variable independent 
from (y4(n))^gpj. Let 0 be topical from to M. Assume one of the following 
conditions: 

i) sup 3 , |0(y4(l)x) — 0 (a:)| has a second moment, 

a) ^(1) 0 has a 4 -|- e-th moment and has a 2 + e-th moment. 

Then there exists > 0 such that weakly to a random 

vector whose coordinates are equal and have law A/'(0,(j^). 


In the first case, or if A{1) 0 has a 6 + e-th moment and X^ has a 3 -|- e-th 
moment, then 

• cr^ = lim (0 {x{n, X°)) — 

• a = 0 iff there is a 9 E Topd with rank 1 such that for any A E Sa and 
any 9' E Ta with rank 1, 9A9' = 99' -\- ■y 1. 

Remark 2.1. According to lemma [3121 if there is such a 9, then every 9 E Ta 
with rank 1 has this property. 

Theorem 2.2 (CLT with rate). Let (A(? 7 ,))^gj^ be an i.i.d sequence of topi¬ 
cal operators with the MLP property and X^ an M.'^-valued random variable 
independent from Let ej) be topical from to M. Assume one of 

the following conditions: 

i) sup 3, |0(A(l)a;) — 0(a:)| has an l-th moment with I > 3, 


a) A(l) 0 has an l-th moment, with I > 6 . 

Ifa^ > 0 in theorem [13 then there exists C > 0 such that for every initial 
condition X^ with an l-th moment, we have 


sup |P[0(x(n, 2^°)) — ny — <p{X^) < auy/n] — A/'(0,1)(] — 00 , u]) | 


< 


C 1 +E 


||X^ 


n 


( 2 ) 
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sup 

ueR'i 


[x{n, X°) — ri'jl < au^/n] — M{0,1 )(] — oo, minuj]) 


< 


Cil+E 


ll^°llL 


E 


p(i)oii: 


n2(i+i) 


Definition 2.1. We say that the sequence ( 74 (n))^gj^ is algebraically arith¬ 
metic if there are a, fe G M and a 6 * G T opd with rank 1 such that for any 
A E Sa and any 6' G Ta with rank 1, 


{OAe' - ee'){w^) c (a + 6z)i. 


(3) 


Otherwise the sequence is algebraically non arithmetic. 

Remark 2 . 2 . According to lemmaif there is such a 9, then every 6 E Ta 
with rank 1 has this property. 

Moreover, for any 6, 6' E Topd with rank 1 and any A E Topd, the function 
6A6' — 66' is constant with value in Ml. 


Theorem 2.3 (LLT). Let (A(n))^gj^ be an i.i.d sequence of topical operators 
with the MLR property and an Ef'-valued random variable independent 
from (A(n))^gpj. Let cf be topical from M'^ to M. Assume one of the following 
conditions: 


i) supj, \(j){A{l)x) — (j){x)\ has a second moment, a > 0, and the system is 
non arithmetic 


a) A(l) 0 has a 4 -I- e-th moment, G L°° and the sequence is 

algebraically non arithmetic. 

Then a > 0 and there exists a a-finite measure a on M'^, so that for any 
continuous function h with compact support, we have: 


lim sup 

tiSlR 


crV27rnE \h {x{n, X^) — nyl — ul)] — E 


(u + cl,(xLf 


g 2ncr2 

a{Ji) 


0 . 


Moreover the image of a by the function x i—>■ (x, 0(x)) is the product of the 
invariant probability measure on by the Lebesgue measure. 


Remark 2.3. Like in the usual LLT, this theorem says that the probability 
for x{n,X^) to fall in a box decreases like ■^. To replace the continuous 
functions by indicator functions of the box, we need to know more about the 
invariant probability measure on In particular, numerical simulations 

show that some hyperplanes may have a weight for this probability measure, 
so those hyperplanes could not intersect the boundary of the box. 

















The algebraic non arithmeticity is optimal in the following sense: 


Proposition 2.1. If the conclusion of theorem lEM is true, then 
is algebraically non arithmetic. 

Theorem 2.4 (Renewal theorem). Let be an i.i.d sequence of top¬ 

ical operators with the MLP property and an M.'^-valued random variable 
independent from Assume that there is a topical (f from to 

R such that sup 3 , | 0 ( 74 (l)a;) — 0(a;)| has a second moment. We denote by a 
the same measure as in theorem, fS.A //y > 0 and the system is weakly non 
arithmetic, then for any function h continuous with compact support and any 
initial condition , we have: 

lim \h ix{n,X^) — al)] = 0, 

n>l 

lim y^E [h (^x(n,X^) — al)] = . 

n>l ^ 

Remark 2.4. The vector 1 gives the average direction in which x{n,X^) is 
going to infinity. Like in the usual renewal theorem, this theorem says that 
the average number of x{n, X°) falling in a box is asymptotically proportional 
to the length of this box, when the box is going to infinity in that direction. 
Like in the LLT, to replace the continuous functions by indicator functions 
of the box, we need to know more about the invariant probability measure 

on PRLx- 

Theorem 2.5 (Large deviations). Let (y 4 (? 7 ,))^gj^ be an i.i.d sequence of top¬ 
ical operators with the MLP property and X^ an M.'^-valued random variable 
independent from (y4(n))^gpj. Let cf) be topical from R'^ to R. Assume that 
sup 3 , \(l){A{l)x) — 0(x)| has an exponential moment, and that > 0 in theo¬ 
rem [O Then, there exists a non negative strictly convex function c, defined 
on a neighborhood of 0 and vanishing only at 0 such that for any bounded 
initial condition X^ and any e > 0 small enough we have: 

lim — In (P [0 {x{n, X°)) — ny > ne]) = —c(e), 
lim — In (P [0 ix{n, X°)) — ny < —nel) = —c(—e). 

n Tl ^ ^ ^ ^ 
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2.2 Max-plus case 

When the A{n) are (max, +)operators, it is natural to chose 0(x) = maxjXj. 
In this case we get rniuj max* Aij < (j){Ax)—(j){x) < maxjj Aij, so integrability 
condition can be checked on the last two quantities. 

Theorem 2.6 (CLT). Let be an i.i.d sequence of {max, +)operators 

with the MLP property and an W^-valued random variable independent 
from //max^ A{l)ij and min^ maxj A{l)ij have a second moment, 

then there exists > 0 such that for every initial condition 

(i) weakly to a random variable whose coordinates 

are equals and have law N{t),a'^), 

(a) = lim i f (maxjj A(n) ■ ■ ■ A(l).y)^ dF. 

Theorems 12.21 to 12.51 are specialized in the same way: the conclusion is 
valid if 0 (x) = maXjXj and the moment hypothesis on sup 3 , |0(y4(l)x) — 0 (x)| 
is satished by maxij A{l)ij and mmjmaxiA{l)ij. 

In this case, we also get another condition to avoid degeneracy in the CLT. 
To state it, we recall a few definitions and results about (max, +)matrices: 

Definition 2.2. For any k,l,m & N, the product of two matrices A G 
and B e is the matrix AB E dehned by : 

VI < i < /c, VI < j < m, {AB)ij := max Aip + Bpj. 

i<p<i 

If those matrices have no line of — 00 , then the (max, +)operator defined 
by AB is the composition of those dehned by A and the one dehned by B. 

Definition 2.3. A circuit on a directed graph is a closed path on the graph. 
Let A be a square matrix of size d with entries in Mmax- 

i) The graph of A is the directed weighted graph whose nodes are the 
integers from 1 to d and whose arcs are the {i,j) such that Aij > — 00 . 
The weight on {i,j) is Aij. The graph will be denoted by G{A) and the 
set of its elementary circuits by C{A). 

ii) The average weight of a circuit c = (A, •'' An, in+i) (where A = in+i) is 
aw{A,c) := ^Ei=i 

iii) The (max, +)-spectral radiu^ of A is Pmax(A) := maxcgc(yi) aw{A, c). 
^this quantity is the maximal (max, +)-eigenvalue of A, that is 

Pmax(A) = max{A e Rniax|3F G REx\{(-oo)‘^}, Ay = y + Al}. 

See [Gau92| . 
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Theorem 2.7. Assume the hypothesis of theorem \2.(A with 7 = 0 . Then the 
variance in theorem \2 .61 is 0 if and only if {praax{B)\B G T^} = { 0 }. 

Theorem 3.2 of [Mer04j gives a condition to ensure the memory loss prop¬ 
erty. This condition also ensures that there are two matrices in Sa with two 
distinct spectral radius. This proves the following corollary: 

Corollary 2.1. Let the law of A(l) be a probability measure on the set of 
d X d matrices with finite second moment whose support is not included in 
the union of finitely many affine hyperplanes Then x{n,.) satisfies 

the conclusions of theorem with a > 0 . 

We also give a sufficient condition to ensure the algebraic non arithmetic- 
ity: 

Theorem 2.8. Assume the hypothesis of theorem l^.dl ii) except the alge¬ 
braic non arithmeticity and A{n) are (max,-|-)operators. If (y4(n))^gj^ is 
algebraically arithmetic, then there are a, 6 G M such that 

{pmax(5)|i? G SA, Q{B)strongly connected} C a + hij. 

Together with corollary 12.11 this proves that the hypothesis are generic 
in the following sense: 

Corollary 2.2. If the law of ^(1) is a probability measure on the set of 
d X d matrices with 4 -|- e-th moment whose support is not included in the 
union of enumerably many affine hyperplanes ofMf^'^, then x{n,.) satisfies 
the conclusions of theorem \2.A 

2.3 Comments 


The following table sums up the limit theorems. In each situation we assume 
that the sequence (A(n))^gj^ has the memory loss property. 


Theorems: 

.4(1)0 

Moments of 

maxjj A{l)ij and miuj max* A{l)ij 

Additional 

condition 

CLT 

4 + e 

2 


CLT with rate 

6 -|- e 

3 


LIT 

4 + e 

2 

NA 

Renewal 

~ 

2 

NA 

LDP 

~ 

exp 

G 


NA= non arithmeticity 
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Let us first notice that the results of the second column are optimal in 
the sense that their restriction to d = 1 are exactly the usual theorems for 
sum of i.i.d. real variable (except for LDP). 

The results of the second column are stated for (max, +)operators be¬ 
cause the sup 3, \(l){A{l)x) — 0(x)| is not bounded for general topical opera¬ 
tors. For other subclasses of topical operators, one has to choose cf) such that 
sup 3 , |0(y4(l)a:) — 0(x)| is integrable. For instance 0(x) = minjXj is natural 
for (min, -|-) operators. Actually, it should be possible to derive renewal the¬ 
orem and large deviation principle with the method of [HHOd] but this has 
not been written down. 

The results of the first column require stronger integrability conditions 
but they are also better for two reasons: they are true for any topical op¬ 
erators and the algebraic non arithmeticity does not depend on (j). It is 
expressed without introducing the Markovian operator Q although the sys¬ 
tem is algebraically non arithmetic iff Q has an eigenvector with eigenvalue 
with modulus 1. Moreover for (max,-|-)operators the algebraic non arith¬ 
meticity can be deduced from theorem 12.81 

An important case for (max, -|-)operators is when Aij G M’*' U {—cxo} and 
An > 0 because it modelizes situations where Xi{n ,.) is the date of the n- 
th event of type i and the Aij are delays. In this case the integrability of 
maxij A{l)ij and min^ max* A(l).y is equivalent to the integrability of A(l) 0 . 

We mentioned earlier that J. Resing and ah |RdVH'*~9n] obtained a cen¬ 
tral limit theorem. In a sense our result is weaker because the MLP property 
implies that x{n, .) is uniformly <F-recurrent and aperiodic. But our integra¬ 
bility conditions are much weaker and the MLP property is easier to check. 

F. Toomey’s large deviations principle only requires the uniform bound of 
the projective image that is a very strong integrability condition. It suggests 
that the MLP property should not be necessary. But his formulation of the 
LDP is not equivalent to ours and in the (max, -|-)case it needs the hxed 
structure property, that is P(Ajj(l) = — 00 ) G { 0 , 1 }. 

3 Proofs of the limit theorems 

3.1 Prom iterated topical functions to Markov chains 

In this section, we show that the hypothesis of the theorems on our model 
stated in section [2] imply the hypothesis of the general theorems of |HH01j 
et |HH04j . To apply the results of |HH01j . 
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• the space E will be Topd x with the Borel a-algebra, 

• the transition probability Q will be dehned by 

Q{{A,x),D) = ¥{{A{l),M)eD), 

• the Markov chain will be {A{n),x{n — 1,.)), 

• the function ^ will be dehned by x) = (j){Ax) — 4>{x), where ^ is a 
topical function from to M. 

• (t{A) = supa; x)| = supa; \(j){Ax) — (/>(a:)| < oo a.s.. 

With these dehnitions, Sn = XlILi is equal to (j){x{n, W°)) — 0(X°). 
To apply the results of |HHnij . we still need to dehne the space B. 

Definition 3.1. Let C°° be the space of complex valued bounded continuous 
functions on PM^ax- 

Let j be the function from Topd x PM^ax to iPi^max such that j{A, x) = Ax 
and / the function from M™max to dehned by 


1 ( 0 ) = (j)oj. 

We call B°° the image of by I. 

Since H-Hoo) is a Banach space, / is an injection, and ||0 o I\\^ = 
||</’||oo5 {B°°, ||.||oo) is also a Banach space. 

Definition 3.2. The Fourier kernels denoted by Qt or Q{t) are dehned for 
any t G C by 

Qt{x,dy) = E*^^y^Q{x,dy). 

We say that the system is non arithmetic if Q(.) is continuous from M to the 
space of continuous linear operators on and, for any t E M.*, the 
spectral radius of Q(t) is strictly less than 1. 

We say that the system is weakly non arithmetic if Q{.) is continuous 
from M to /Igcx) and, for all t G M*, Id — Q{t), where Id is the identity on 
B°°, is invertible. 

Proposition 3.1. If a (^(1)) has an m-th moment and if A{n) has the MLP 
property, then {Q,f,B°°) satisfies condition H{m) of Moreover, the 

interval Iq in condition {HS) is the whole M and s{Q,B°°) = 1. 

To prove (Lf2), we will use theorem 11.31 In the sequel, z/q will be the law 
of Y in that theorem, that is the (unique) invariant probability measure. 
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Proof of proposition \3.1[ Condition {HI) is trivial, because of the choice 
of 


To check condition {H2), we take z/ := /r 0 z/q. It is Q-invariant by 
dehnition of z/q. This proves (i). To prove {ii) and {in), we investigate the 
iterates of Q. For any 0 e and x G we have: 

\Q^{(j)oj){A,x)-iy{(j)oj)\ = \Q^{(j)){Ax) - iyo{(j))\ 

= j (j) {x{n, Ax)) dF — J (p{Yn)dF 

< ||0||oo2P(3xo,10 7^ :r(n,Xo)), 

If we denote 0 i—^ ^^(0) by N, we obtain: 

\\Q^-N\\<2F{3x,Y,,^x{n,x))^0. (4) 

This proves that the spectral radius r{Q\KerN) is strictly less than 1 and 
that sup„ IIQ^II < oo- Since Q\imN is the identity, dim{ImN) = 1 and 
= KerN © ImN, Q is quasi-compact, so {ii) is checked. Moreover 
s{Q,B°°) = 1. It also proves Ker{\ — Q) C. ImN, which implies {Hi). 

To prove {H3) we set := P^^'^'^\ii{y))^Q{x,dy) and 

(5) 

To prove that is the derivative of Q^, it remains to bound by 

a quantity that tends to zero with h. 

To this aim, we introduce the following function: 

f / : M ^ C 
1 t I—> e** — 1 — if. 


The calculus will be based on the following estimations on /: \f{t)\ < 2t, 
and 1 / 0)1 < 0 . 

Now everything follows from 

Ai^\<l)oj){A,x)= j 0 (:^) {ii{B,Ai)f f {haB,7^))dy{B). 

( 6 ) 

First it implies that 


A^^\(t)oj) <1101 


a\B)\\f{hi{B,.))\\^dy{B). 


(7) 
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Since |/(t)| < 

i^<T(B)‘||/(/!«B,.))ll„ < ^ 0. 

Since |/(t)| < 2t, 

||/(Ae(B,.))IL<2(r‘+'(B). 

When k < m, is integrable so the dominate convergence theorem 
and the last two equations show that 

j <r(Bf\\f(h((B,.))\\^MB) = o(h). (8) 

Finally for any k < m, tends to zero, so is the derivative 

of in t. 

To prove that is continuous, we notice that 

(qS - Qt’) (<t'ojKA,x) = U (BS) g(hHB,M)) dti(B). 

(9) 

where g{t) = e®* — 1. Then we apply the same method as before, replacing 
the estimates on / by \g(t)\ <t to prove the convergence, and by \g(t)\ < 2 
to prove the domination. 

This proves {H3) and the additional assumption of proposition 13.11 □ 

In their article [HHOd] H. Hennion and L. Herve have proved limit theo¬ 
rems for sequences ^(W, Z„_i), where (W)neN is an i.i.d. sequence of Lips- 
chitz operators on a metric space Ai, and is defined by Zn+i = As 

explained in section [T73l we take Ai = W = A{n) and again .^(A, x) = 

(p{Ax) — (p{x). In this case Zn = x{n,X^) and Sn = (j) {x{n, X^)) — 0(X°). 

Moreover in our situation, the W, which are the projective function defined 
by A{n), are 1-Lipschitz. Following the same proof as [HH04j with this ad¬ 
ditional condition, we get the CLT (resp. CLT with rate, LLT) for Sn under 
the hypothesis of theorem 12.II (resp. 12.2112.31) on A(1)0. 

The integrability conditions are weaker than in [HHOd] , because the Lips- 
chitz coefficient is uniformly bounded. The only difference in the proof is the 
Holder inequality of the 4th part of proposition 7.3 of |HH04j : the exponents 
in the inequality should be changed to 1 and cx). 

Let us give the notations of |HHn4j we need to state the results. G is the 
semi-group of the operators on Ai, 1-Lipschitz for distance 6. For a fixed 
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Xq G A4, every rj > 1 and every n G N, we set Airj = ^^o)] and 

Cn = ^[c{Yn ■ ■ ■ Vi)], where c(.) is the Lipschitz coefficient. 

When there is an iV G N such that Cn < 1, there is a Aq g] 0, 1[, such 
that jQc{g) (1 + XoS{gxo, xo))'^'^ dfi*^ {g) < 1. We chose one such Aq and set 
the following notations: 


(i) Br^ is the set of functions / from W1 to C such that rnr,{f) < oo, with 
norm ||/||^ = |/|^ + where 


I/I. 

m^if) 


sup 

X 

sup 


(1 + Ao(5(x,Xo))^+^’ 

_ \f{x)- f{y)\ _ 

6{x,y) (1 + Ao(5(a:,a;o))’' (1 + Xo6{y,xo)y' 


(ii) We say that the system is ? 7 -non arithmetic if there is no t G M\{0}, no 
p G C, and no tc G with non-zero constant modulus on the support 
SuQ of the invariant probability measure uq such that |p| = 1 and for all 
n G N, we have 

e^*^-w{Zn) = p^w{Zo)¥ - a.s., ( 10 ) 

when Zq has law z/q. 

Remark 3.1 (non arithmeticity). In the first frame the non arithmeticity 
condition is about the spectral radius of Qt- Here we work with the associated 
Pt that acts on M. instead oi G x M. (cf. |HH04j h If Pt is quasi-compact, 
then the spectral radius r{Pt) is 1 iff Pt has an eigenvalue p with modulus 1. 
It is shown in proposition 9.1’ of [HHOd] that if r{Pt) = 1, then Pt is quasi¬ 
compact as an operator on Bjj and that an eigenvector w with eigenvalue p 
satisfies equation CB. 

Proposition 3.2. 

1. If A{1)0 has an rj-th moment, with rj G then Air] < oo. If the 
sequence has the MLP property, then there is Uq E N such that Cn^ < 1. 
If has an rj-th moment y G M"*", then f i—^ E[/(X°)] is continuous 
on Brj. 

2. Algebraic non arithmeticity implies rj-non arithmeticity for any rj > 0. 

The first part of the proposition is obvious. The second part relies on the 
next two lemma that will be proved after the proposition; 

Lemma 3.1. The support of the invariant measure z/q is 
Si^Q := {91\9 G Ta, 9 with rank 1}. 
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Lemma 3.2. If equation ^ is satisfied by some 6 with rank 1, any A E Sa 
and any 9' E Ta with rank 1, it is satisfied by any 6 eTa with rank 1. 

Proof of proposition 13.21 Let us assume that the system is r^-arithmetic. Then 
there are w E and t,a eM such that for /r-almost every A and uq almost 
every x, we have: 

( 11 ) 

Since all functions in this equation are continuous, it is true for x E and 
A eTa- Since is Ta invariant, we iterate equation flTTl) and get 

^imTx)-^{x))^(7^^ = ( 12 ) 

where T E Ta and ht is the number of operators of Sa to be composed to 
obtain T. 

Because of the MLP property, there is a. 9 E Ta with rank 1. For any 
A E Sa, 9A E Ta, so we apply equation flT^ for T = 9A and T = 9 and 
divide the first equation by the second one. Since ngA = ng + 1 and 9Ax = 9x 
, we get 

^it{(l){6Ax)—(p{9x)) _ ^ita 

Setting 6 = ^, it means that fiifiAx) — (j){9x) E a + bT,. Since 9 has 
rank one, {9Ax — 9x) E Ml, so 9Ax — 9x E (a + 6Z)1, and the algebraic 
arithmeticity follows by lemma 13.11 □ 

Proof of lemma 1X71 By theorem 11.31 there is sequence of random variables 
Yn with law z/q, such that Yn = A{n) ■ ■ ■ A{1)Y. Let K be a compact subset 
of Mmax such that Y E K with positive probability. 

For any 9 E Ta and any e > 0, the set V of topical functions A such that 
6{Ax, 9x) < e for all T G iF is a neighborhood of 9. Therefore the probability 
for A{ng) ■ ■ ■ A(l) to be in V is positive and by independence of F, we have: 

P[F e K,A{ng)---A{l) G F] > 0. 

Since 91 = 9Y, this means that with positive probability, 

6{Yn„Wl) =6{A{ne)---A{l)Y,W) < e. 


so 91 E Syg. 

This proves that {91\9 E Ta, 9 with rank 1} C Syg. 

In |Mai97j . is obtained as the law oi Z = lim„ ^4(1) ■ ■ ■ y4(n)l. In¬ 
deed, the MLP property and the Poincare recurrence theorem ensure that 
there are almost surely M and N such that A[N) ■ ■ ■ A{N + M) has rank 1. 
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Therefore, for n > iV + M, A(l) ■ ■ ■ A{n)l = A(l) ■ ■ ■ A{N + M)1 = Z . But 
A(l) ■ ■ ■ A[N + M) G Ta almost surely, so Z G {6*1|6* G Ta-,0 with rank 1} 
almost surely and C {91\9 G Ta,9 with rank 1}. □ 

Proof of lemma [Ql We assume that equation ([3]) is satished by 6 * = 6 *i, any 
A E Sa and any 9' G Ta with rank 1. 

Let Ai, ■ ■ ■, An E Sa, such that 92 = Ai ■ ■ ■ An has rank 1. For any i < n, 
Ai - ■ ■ An9' has rank 1, so (6*1^4* • ■ ■ An9' — A+i ■ ■ ■ C (a + 6Z)l. 

Summing these inclusions for z = 1 to z = n, we get {9i929' — 9i9'){W^) C 
(na + 6 Z )1 and 

{{9i 92A9' - 9iA9') - [9x929' - 9x9')) (M'^) C hlA. (13) 

Now we write 929' as 

Q^Q' = 9x9' + [9x929' - 9x9') - [9x929' - 929'). 

The last part does not depend on 9', so replacing 9' by A9' and subtracting 
the hrst version, we get: 

92A9' - 929' = 9xA9' - 9x9' + [[9x92A9' - 9xA9') - [9x929' - 9x9')). 
With equation flT^ . this proves equation ([3]) for 6 * = 6 * 2 . 

□ 

3.2 Prom Markov chains to iterated topical functions 

Propositions 13.11 and 13.21 prove that under the hypothesis of section [2] the 
conclusions of the theorems of [HHOl] and [HHOd] are true. This gives results 
about the convergence of (0 [x[n,X^)) — (j)[X^) — nh'[^),x[n, X^)). 

When X^ has law uq, the sequence (^(n),^(n,is stationary, so 
it follows from Birkhoff theorem that 7 = / f,[A,x)dh'o[x)dfi[A) = i'[f,). 

The following lemma will be useful to go back to a;(? 7 ,,.). 

Lemma 3.3. If (p is a topical function from to M, the function ip ■. x 
[(p[x),x) is a Lipschitz homeomorphism with Lipschitz inverse from onto 

R X 

Proof. Let [t,x) be an element of R x Then ip[y) = [t,x) if and only 

if there is an a G R such that y = x + al and (p[x) + a = t. So the equation 
has exactly one solution y = x + [t — (p[x))l and ip is invertible. 

It is well known that topical functions are Lipschitz, and the projection 
is linear, so it is Lipschitz and so \s ip. 
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For any x,y ^ we have x < y + maxj(a:i — yi)l, so (j){x) — (j){y) < 
ma.Xi{xi — yi). Therefore, for any 1 < z < h, we have 

(t){x) - (t){y) - {xi - yi) < max(a:i - yt) - min(a:j - yi) = 6{x,y). 

i i 

Permuting x an y, we see that: 


|0(x) -0(z/) - {xi-yi)\ < S{x,y). 

Therefore \xi — yi\ < \4>{x) — </>(|/)| + S(x,y) and ^jJ~^ is Lipschitz. 


(14) 

□ 


Proof of theorem \2. 1[ Without lost of generality, we assume that 7 = 0 . The¬ 
orem A of [HHOl] and proposition 13.11 or theorem A of |HH04j and proposi¬ 
tion prove that ) converges to A/'(0, a^), which means that 

<l>{x{n,xj^- 4 >{x ) ^ converges to the limit specihed in theorem 12. II . We just es¬ 
timate the difference between the converging sequence and the one we want 
to converge: 


An : = 


x{n, A°) 
Vn 


(f{x{n,X^))-(f{X^) 


|0(X°)| ^ \x{n,X^)\r 
Vn ^/n 


(15) 

Each term of the last sum is a weakly converging sequence divided by y/n 
so it converges to zero in probability. This proves that An converges to zero 
in probability, which ensures the convergence of ^ to the Gaussian law. 


The expression of is the direct consequence of theorems A of |HHnij 
or theorem S of [HH04] . 


If cr = 0, then again by theorem A of |HH01j or S of [HH04j . there is a 
continuous function on PM^ax such that 

(t){Ax) - (t){x) = f{x) - f{Ax) (16) 

for yU-almost every A and i^o-almost every x. Since all functions are continuous 
in this equation, (lT 6 l) is true for every A ^ Sa and x G By induction 

we get it for A E Ta and ii 9 E Ta has rank 1 and x E 9Ax = 9x, so 

(j){9Ax) = (j){9x). 

Since 9Ax — 9x E Ml, this means that 9Ax = 9x. By lemma IHTTl it proves 
that 9A9' = 99' for any 9^ 9' E Ta with rank 1 and A E Sa- 


Conversely, let us assume there is 9 with rank one such that for any 
9' E Ta with rank 1 and A E Sa, 'we have: 

9A9' = 99'. (17) 
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By lemma [H721 applied with a = fe = 0, it is true for any 6,6' & Ta with rank 1, 
and any A E Sa and by induction, equation ffT7|) is still true for A E Ta- 
Therefore, for any m eN and n > m + 1 and any 6' E Ta with rank 1, if 
A{n) ■ ■ ■ A{n — m + 1) has rank 1 ,then x{n, 6'1) = A{n) ■ ■ ■ A{n — m + 1)6'1 
and for any N eN 

F{\\x{n,6'l)\\oo<N) 

> P iA{n) ■ ■ ■ A{n — m + l)has rank 1, || A(n) ■ ■ ■ A(n — m + l)6*'l||oo < N) 

> P (^(m) ■ ■ ■ 74(l)has rank 1, ||74(m) ■ ■ ■ A(l)6*'l||oo <-/V). (18) 


We £x a 6 ^' G Ta with rank one. The MLP property says there is an m 
such that ¥{A{m) ■ ■ ■ 74 (l)has rank 1) > 0. Therefore, there is an iV G N 
such that the right member of (lT 8 l) is a positive number we denote by /3. 

Equation (lT 8 l) now implies that for any e > 0, if n > max(m, 
then P(||;^x(n, 6*'l)||oo < e) > /3, so A/’(0, (T^)[—e, e] > /3. When e tends to 
zero, we get that A/’(0, cr^)({0}) > P > 0, which is true only if a = 0. □ 

Proof of theorem \2.A Without loss of generality, we assume that 7 = 0. 
Equation ((21) follows from theorem B of [HHOl] and proposition 13.11 or from 
theorem B of [HHOdj and proposition 13.21 

The only fact to check is that the initial condition dehnes a continuous 
linear form on Bn, with norm at most C (l + E(||X°||(,^), that is for any 
f E Bn, 'we have: 

|E(/(A-»))|<C(1+E(||A'»|7)) ll/ll,. 

It easily follows from the fact that |/(a;)| < ||/||,,(1 + and l + rj <1. 

Taking ?/ = 0 in (|T4l) . we get |0(x) — Xi\ < \x\'p. Together with ((T5l) it 
proves that for any u G and any e > 0 


[x{n,X^) < auy/n] 
< P 


< P 


minxiin, X^) < amin.Ui^/n 
i i 

4>{x{n,X^)) < (a minwj + 2e)A/n 


+ P 


| 0 (XO)| 


> e 


n 


+ P 


\x{n,X^)\v 


> e 


n 


,,, 2e,, C E(|a;(n,X0)|U 

< Ar(0,1)( - cx),mmMi H-) + —H- - — 

^ a P\/ny 


+ 


(e 0 r)' 


< A/'(0,1)(] - cx),minMj) + ^ H-h , 

i ^/n a PV'^y 


C 2e E(|T(n,X0)|),) 




(19) 
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Conversely, 

\x{n, X°) < auy/ri] 
> 


> P 


0 (a;(n, X°)) < a minui^/n 
i 

0 (x(n, X°)) < (a minMj — 2 e) 


P 


| 0 (X°)| 


> e 


n 


P 


|x(n,X°)|p 


> e 


n 


> ■ 2^n C e(W^°)I') E(|$(n,.Y»)|y 

> A/(0,1)(J — cxD, minMj-J)- -= — 

^ (T xj Ti 


> A/”(0,1)(] — cxD, mintij])— 


C 


n 


2 e 

a 


[ey/ny 

E(\cP{X^)f 


{ey/^y 


{ey/nY 

E(|x(n,X°)|^) 

{ey/nY 


( 20 ) 


Taking e = n~ in flT^ and (12(11) will conclnde the proof of theo¬ 
rem [212] if we can show that E (|x(n, X°)|p) is bonnded nniformly in n and 
X^. Withont loss of generality, we assnme X^ = 0. 

For no G M, we take a > (P [A(no) ■ • ■ A(l) has not rank 1 Bnt if 

A{n) ■ ■ ■ A{m) has not rank 1, then for any integer less than , the oper¬ 

ator A(l + mo) ■ ■ ■ A((i-|-l)no) has not rank 1 either. From the independence 
of the A(n), we dednce 


n—m—no 


P {A{n) ■ ■ ■ A{m + 1) has not rank 1 ) < a 

We estimate 5 {A{n) ■ ■ ■ A{m + 1)0, A{n) ■ ■ ■ A{nQ + 1 + m)0): it is 0 when 
yl(n-|-m) ■ ■ ■ y4(no+l+m) has rank 1 , and it is always less than <5 (A(no -l- rn) ■ ■ ■ A{m + 1 ) 0 , 0 ), 
that is less than p^s not rank 1 } + m) ■ ■ ■ kl(m + 1 ) 0 |^, 

where 11 denotes the indicator fnnction. Therefore, we have for any n > m+UQ 

E [5* {A{n) ■ ■ ■ A{m)0, A{n) ■ ■ ■ A(no -|- 1 -|- m)0)] 


< E 


1 


= a 


{A{n)-A{no+m+i) has not rank 1 } l^(^o + m) ■ ■ ■ A{m + l)0|p 

-""°E[|kl(no)---kl(l)0|;,l . (21) 


|A(no) ■ ■ ■ A(l)0|p 
Let n = qno -|- r be the Enclidean division of n by riQ. Then we have 

|a;(n, 0)|p = 5 (A(n) ■ ■ ■ A(1)0, 0) 

<? 

< 6 {A{n) ■ ■ ■ A{ino + 1)0, A{n) ■ ■ ■ A{{i — l)no -|- 1)0) 

2=1 


-|-(5 {A{n) ■ ■ ■ A{n — r + 1 ) 0 , 0 ). 
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Therefore we have: 




- [l^K)---^(l)0| 

i=l 


i/i 


E 




in 


(22) 


We apply this decomposition again (with n = r, no = 1 and a = 1), to check 
that 


E 


|kl(r)---kl(l)0| 


V 


in 


< r (E 


|kl(l)0| 


V 


in 


< no (E 


|kl(l)0| 


V 


in 


It follows from the MLP property, that there is no G N such that a < 1. 
Introducing the last equation in equation fl^ . we see that 


E 


\x{n, 0)1 


V 


i/i 


< 1 + 


a 


— ‘ItlqI 


1 — 


no 


E 


|kl(l)0| 


V 


in 


□ 


To go from the abstract LLT and renewal theorem to ours, we will use 
the following classical approximation lemma. 

Lemma 3.4. Let h be a continuous function with compact support from 
X Then there are two continuous functions fo and go with compact 
support in and Mf respectively, so that for any e > 0, there are fi and Pi 
continuous functions with compact support satisfying: 

Wx e M“,?/ e \h{x,y) -Y^fi{x)gi{y)\ < efo{x)go{y). 

i 


In the sequel, we denote by C the Lebesgue measure. 

Proof theorem 1^.31 By theorem 12.11 the algebraic non arithmeticity ensures 
that cr > 0. 

We apply proposition 13.II and theorem C of |HH01j or proposition 13.21 and 
theorem C of |HHn4j . This proves that, if G Cc(lR) and if / is a bounded 
Lipschitz function on then for any x^ G M'^, 


lim sup 

ueR 


aVrimE (/ (a;(n, a;°)) g (</> (x(n, x^)) — (j) (a;°) — u)) — e 


Moreover these convergences are uniform in x^, because 6-^ 
independently of xq as a linear form on B°° and is in a disk of Bri 


(23) 
is bounded 
with center 
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0 and radius ||X°||oo if \x^\v < 11I loo- The uniformity allows us to take any 
random initial condition and get 


lim sup 

" ueK 


aV^nnE [f (x(n, X°)) g (0 (x(n, — «)] — E 


( u +</,( X °))2 
g 2710-2 


(24) 


But the density of bounded Lipschitz functions in ||.||oo) al¬ 

lows us to take / and g continuous functions with compact support in equa¬ 
tion fl24p . Now, it follows from lemma [H~4l that for any h continuous with 
compact support 


0 . 


lim sup 

ueR 


aV^nnE [h (a;(n, 0 (x(n, 


— «)] — E 


(u+^(x0)p 


g 2710-2 

oo (8) £(h) 


(25) 

According to lemma 13.31 the function $ : a; i—*• (0, x) is Lipschitz with a 
Lipschitz inverse, therefore h has compact support iff h o <h does. Since 
a; -I- m1 = x, this concludes the proof. □ 


0 . 


Proof of proposition \2.1[ Assume the sequence of random variables is alge¬ 
braically arithmetic and the conclusion of theorem 12.31 holds. 

There are a, 6 G M and 6 with rank 1, such that every A ^ Sa and 6' G Ta 
with rank 1 satisfy equation ([3]). We set t ^ if 6 7 ^ 0 and t = 1 otherwise, 
and for any x G EA, w{x) = _ Equation ([3]) implies that, for any 

A G Sa, y G M'^, and any O' G Ta with rank 1: 




We chose y such that <f>{6'y) = 0. By induction, we get 

(26) 

For any / : M 1 —^ M and g : PE^ax ® continuous with compact supports, 
the conclusion of theorem l2.3l for h dehned by h{x) = f{(j){x)){gw){x) is that: 

aV^E [/(0(x(n, 9'y)))g(x{n, 9'y))w(x{n, 9'y))] C{f)uo{gw). 

Together with equation (|26|) . it means 

P^''''w{e'y)aV^E [e"**-/(0(x(n, e'y)))g{x{n, O'y))] C{f)vo{gw) (27) 

But conclusion of theorem 12.31 for h dehned by h{x) = (/e“**')(0(x))5f(T) 
is that: 


(T\/2TmE[e **-/(0(x(n, 0'|/)))^(x(n, 0'?/))] ^ C{fe '^')vo{g) (28) 
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Equations IHT} and fl^ together imply that ta E 2tt'L and that 
w{9'y)C{fe"^-)i'Q{g) = C{f)vQ{gw). 

The right side of the equation does not depend on 9' so by lemma 13.11 w is 
constant on this proves I'a^gw) = w{9'y), so £(/e“®*') = C{f) that is 
e®*' = 1 or t = 0. This is a contradiction, which concludes the proof. □ 

Proof of theorem \2.4\ Applying proposition 13.II and theorem D of |HHnij . we 
have that, ii g E Cc(M) and if / is a bounded Lipschitz function on 
then for any G M'^, 

lim ^ E [/ (0 (a;(n, x°)) — 0 (a;°) — a) g [x{n, a:°))] = 0, 

n>l 

lim ^ E [/ (0 [x{n, x^)) — (j) (a;°) — a) g [x{n, a:°))] = ^ 

oo y 

Moreover these convergences are uniform in x°, because 5^ is bounded as a 
linear form on . The uniformity allows us to remove the 0 (x°) in the last 
equations and take any random initial condition. The result follows by the 
same successive approximations as in the proof of the LLT. □ 

Proof of theorem \2.5[ Without lost of generality, we can assume that 7 = 0. 

The exponential moment of cr{A) means that there is a 6* > 0 such that 
j < oo_ An easy bound of the norm of dy) inspired by 

the proof of proposition 13.11 ensures that z i—>■ is analytic on the open ball 
with center 0 and radius 9. To prove that it is continuous on the domain 
{\lZz\ < 6^/2}, we apply the same method. 

Now theorem E of |HH01j gives 

lim — InP [0 (a;(n, X°)) — 0(W°) > nel = —c(e). 

Let 0 < 7 < e. For any n > ||0(X°)/?7||oo, we have 

P [0 (x(n, X°)) — 0(X°) > ne] > P [0 {x{n, X°)) > n(e + rj)] , 
which implies that 

lim inf — InP [0 (x{n, — ny > nel > —c(e + rj). 

The same method gives 

lim sup — InP [0 (a;(n, — ny > ne] < —c(e — y). 

n ^ 

By continuity of c, the hrst equality is proved. The second one follows from 
the same method applied to —0 instead of 0. □ 
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3.3 Max-plus case 

Before proving the statements, we recall a few needed definitions and results 
about powers of matrices in the (max, +)algebra. 

Definition 3.3. 1. The critical graph of A is obtained from Q{A) by 

keeping only nodes and arcs belonging to circuits with average weight 
Pmax(^)- If will be denoted by G^{A). 

2. The cyclicity of a graph is the greatest common divisor of the length 
of its circuits if it is strongly connected (that is if any node can be 
reached from any other). Otherwise it is the least common multiple of 
the cyclicities of its strongly connected components. The cyclicity of A 
is that of G^{A) and is denoted by c{A). 

Remark 3.2. Interpretation of powers with G{A). 

If ■ ■ An) is a path on G{A), its weight is so that 

iA"’).j is the maximum of the weights of length n paths from i to j. 

Theorem 3.1 ( |(JDQV83| ). Assume Q(A) is strongly connected, Pmax(^) = 0. 
Then the sequence is ultimately periodic and the ultimate period is 

the cyclicity of A. 

Proof of theorem \2^ Suppose that a = 0. By proposition 13.II we may apply 
theorem A of [HHOlj . The third point of the theorem says that there exists 
a bounded Lipschitz function / such that for //-almost every {B,x): 

ma.x{Bx)i — maxa:, = f(x) — f{Bx) (29) 

i i 

Since all functions in that equation are continuous, every B E Sa and 
X G Supp{h'o) satisfy equation fl29|) . li B E Sa and x E Suppi^uo), then 
Bx E Suppivo), so by induction equation fl^ is satished by any B in Ta- 
Since for B E Ta, B'^ E Ta, maxj B^Xi is bounded. But there exists a k 
such that c{B)prni,AB) = Bf^\ so max* > nc{B)p^^AB) + Xk 

and pmax(5) < 0. 

Since every path on Q [B) can be split into a path with length at most d 
and closed paths whose average length are at most Pmax(-B), we have: 

max(i?"x)i < (n — d)p^a,x{B) + d max \Bij \ + max a:,, 

i Bij>—oo i 

therefore Pmax(-B) > 0. 

So (T = 0 implies that Vi? E Ta, Pm&AB) = 0. 
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Conversely, if Pmax(-B) = 0 for every B E Ta, then 
maxa;j(?7,, 0) = max (^(n) ■ ■ ■ A(l)).. 

i ij ^ 

> Pmax (-4(n) ■ ■ ■ A(l)) = 0 a.s.. 

Therefore A/'(0, (j^)(M+) > 1, and a = 0. □ 

Proof of theorem We assnme the system is algebraically arithmetic. Then 
there are a, 6 G M and 9 E Ta snch that for any A E Sa and 6' E Ta with 
rank 1, we have: (9A9' — 99')(R'^) C (a + bZ)l. Replacing 9' by AP9', we get 
{9A'^^^9' — 9A'^9'){W^) C (a + bZ)l and by indnction 

{9A^+’^9' - 9A^9')(R^) c {ka + bZ)l (30) 

From now on, we assnme that Q{A) is strongly connected. The matrix A 
dehned by Aij = Aij — pmax(^) satisfy pmax = 0 and has a strongly 
connected graph. Therefore, by theorem 13.11 there is an n snch that for any 
indices i,j, = R”-. 

Since for any n E N, Afj = Afj + npniax(^), it means that = 

Aff'^^ + c{A)prai,^{A), and {9A^+^^^^9')ij-{9A^9')ij = c(A)pmax(^)- Together 
with eqnation (l30l) . it says that Pmax(^) E a + C a + ^Z, which 

conclndes the proof. □ 
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